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ABSTRACT 


Complementary binary code sequences were invented by 
M. J. E. Golay in the investigation of infrared multislit 
spectrometry. Complementary coding sequences have the prop¬ 
erty of an infinite correlation peak to peak ambiguity ratio 
when detected with a matched filter. 

Cooperative or totally orthogonal complementary code 
pairs are two sets of complementary pairs such that the cross 
correlation is zero in every position. A proof ’is given that 
every complementary pair has two totally orthogonal pairs, 
i.e., one the complement of the other. A proof that these 
pairs are the only pairs is also given. 

A communication system involving complementary code 
binary sequences is simulated on the hybrid computer and com¬ 
pared with an ideal receiver for an uncoded signal. By using 
both the totally orthogonal code and time shifting, a method 
of horizontal multiplexing of binary coded information is 
proposed and evaluated. Various other complementary coding 
systems and possible uses are discussed. 
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I. 


INTRODUCTION 


As defined by M. J. E. Golay [Ref. 1], a set of comple¬ 
mentary sequences is a pair of equally long, finite sequences 
of two kinds of elements which have the property that the 
number of like elements with any given separation in one se¬ 
ries is equal to the number of pairs of unlike elements with 
the same separation in the other series for all separations. 
Example: 

A. 00010010 

B. 00011101 

Consider series A at a separation of one element, count¬ 
ing the number of likes: 

00010010 
lluuluu 3 likes 

Next, counting the nximber of unlikes in series B with a sepa¬ 
ration of one element: 

00011101 
llulluu 3 unlikes 

Similarly for longer separations the number of likes in A and 
the number of unlikes in B are listed below: 

Separation Distance Number of Likes A Number of Unlikes B 

2 3 3 

3 4 4 

4 2 2 

5 2 2 

6 11 

7 11 
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The basic property is also expressed in autocorrelative 
terms. Let the a^ and elements (i=l,2,•••,n) of two n- 
long complementary sequences be either +1 or -1, and let the 
respective autocorrelative series be described by the series 
Cj and dj respectively, with subscript ranging over the 
integers from -n+1 to n-1, and defined by the respective 
equations. 

for j ^ 0 


for j <0 

for j 7 ^ 0 

The first use of these complementary sequences was by 
Golay in the field of multislit spectrometry. Since then, 
they have been proposed for radar and sonar. Golay mentioned 
in 1961 that these sequences might be applied to communica¬ 
tion systems [Ref. 1]. 

One problem in radar and sonar systems is to get a maxi¬ 
mum amount of energy radiated for long-range detection, but 
at the same time keep good range resolution. Three basic 
methods are available; FM modulated wave, integration of 
successive pulses,’and intrapulse coding. Of these only 
intrapulse coding will be discussed. 


i=n -3 

'^3 “ ill 


i=n-: 

d. = y b.b.^. 

D i£l ^ ^+3 


Cj = c_^ 

d. = d . 
3 ~3 


Then: c. + d. = 0 

c + d = 2n 
o o 
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Figure 1. Matched filter detection of a random code. 



Figure 2. 


Matched filter detection 


of an 


optimal noise 


code. 
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Using intrapulse coding with a matched filter is illus¬ 


trated in Fig. 1. Note that there is a certain amount of 
coded"signal output from the filter which cannot be complete¬ 
ly eliminated. However, by proper choice of the coded pulse 
with its matched filter it can be minimized. See Fig. 2 for 
an illustration of an "ideal" code filter. It is possible to 
use two complementary codes to eliminate clutter completely. 
An example of this is illustrated in Fig. 3. 

The basic way to make a complementary pair communication 
system is to have two symbols representing a 1 and -1, and 
send either one code or the complement of the code and to re¬ 
ceive the signal in a matched filter. See Fig. 4. 

Example; 


Signal for +1 
Ch 1 11 1-1 

Ch 2 1 1-1 1 


Filter 
1 1 1-1 
1 1-1 1 


Signal for -1 
Ch 1 -l-l-l+l 

Ch 2 -l-l+l-l 


Filter 
1 1 1-1 
1 1-1 1 


Correlation 
-10141 0-1 
1 0-1 4-1 0 1 
0 0 0 8 0 0 0 
Signal decoded as +1 

1 0-1-4-1 0 1 
-1 0 1-4 1 0-1 
0 0 0-8 000 


Signal decoded as -1 

The property that eliminates clutter from the radar sys¬ 
tem and allows the improvement in resolution from a series of 
pulses is the orthogonality of the code with a time shift of 
itself. In a communication system this would allow the 
sender to start the next message one bit behind the previous 
one and by superposition of the linear system it can be seen 
that these bits will be separated at the receiver end. In the 
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Figure 3. Matched Filter Detection of a Complementary 
Code of Length 8. 
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example below a message is sent with the first bit on the 
right, and in the correlation the first bit is on the right; 
Example; 

Encoding Channel 1 Code 111-1 

Message 

- 111 - 1-1 

-l-l-l+l first bit 
-l-l-l+l second 
+1+1+1—1 etc. 


+1+1+1-1 Channel 1 filter Correlation 

-l-l-l+l 


-1 0+1+2-2-3 0+1 

+1+1+1—1 

+1-1-2-2+5+4-3-6-2+1+1 

Encoding Channel 2 

Code 11-11 

filter 


Message 

+1+1-1+1 

-1+1+2-6+3+4-5-2+2-1-1 

-111-1-1 


0 0 0-8+8+8-8-8 000 

-l-l+l-l 

-l-l+l-l 


Message Out 

+1+1-1+1 

+1+1-1+1 

-l-l+l-l 

-1 0+3-2-2+1 0-1 

This could be thought of as five messages added together: 

First 

Message 

000 0-1 

Second 

Message 

0 0 0-1 0 

Third 

Message 

0 0+1 0 0 

Fourth 

Message 

0+1 000 

Fifth 

Message 

-10000 


Sum 

-l+l+l-l-l 

Using superposition 

the s ame 

results are obtained 
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From message 1 
2 

3 

4 


000000 0-8 000 
00000 0-8 0000 
0 0 0 0 0+8 0 0 0 0 0 
000 0+8 000000 


5 00 0-8 0000000 

0 0 0 - 8 + 8 + 8 - 8-8 000 

Two complementary codes are defined to be totally orthog¬ 


onal if their cross correlation is zero in every position. 
Although it was previously found in Ref. 3 that these total¬ 
ly orthogonal codes could be generated from shorter codes for 
for lengths which were a power of two, in this thesis the 
author will show that for every complementary code there are 
two codes that are totally orthogonal to the original code. 

An example of two totally orthogonal codes is shown below. 

Example: 


Original 

Code 


Totally 

Orthogonal 

Code 


Cross Correlation 


Channel 1 +1+1+1-1 -l+l-l-l 


- 1 - 2-1 0 + 1 - 2+1 


Channel 2 +1+1-1+1 -l+l+l+l 


+ 1 + 2+1 0 - 1 + 2-1 


0 0 0 0 0 0 0 


Since the totally orthogonal code is also a complementary 
pair, messages can be encoded both on the original code and 
on the orthogonal code and its complement. The results can 
be added algebraically and transmitted. Thus the information 
on a pair of bandwidths is doubled (See Fig. 5). 

By using the results of sending messages delayed only one 
bit and by also using the orthogonal complementary pair, it 
is possible to send two signals per bit rate on the two 
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channels. If the code length is n, then each bit on each 
channel would have components from 2n messages. 

A detailed simulation of the above horizontal multiplex¬ 
ing was done on the XDS-9300 and CI-5000 hybrid computer, and 
comparison made with an ideal uncoded system. A number of 
variations of communication systems are proposed and dis¬ 
cussed based on complementary binary sequences. 
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? codes 
Codes 


Of 



different lengths is principally one of finding kernels of 
the shortest length that can be confined into the desired 
length. Kernels have been found for lengths 2, 10, and 26 
by Golay [Ref. IJ . Jauregui [Ref. 2] verified by an exhaus¬ 
tive computer search that the kernel found for length 26 and 
its transformations were the only codes to exist at that 

length. Summarised below is the current status of the search 
for kernels. 

Code Length Nir^er of Kernels 

transformations) 


None exist [Ref. 1] 

26 ^ 

None found after extensive but not 
exhaustive search [Ref. 2] . 

It has been hypothesized that none 
exist [Ref. 2] . 

B. TOTALLY ORTHOGONAL COMPLEMENTARY SEQUENCES 

“Two pairs of complementary sequences A^, A , B are 

defined to be "totally orthogonal" if a^^ correlated with 
plus B^ correlated with B^ is zero. More simply stated two 
complementary pair sequences are totally orthogonal if the 
complementary cross correlation is the null vector. Speiser 
and Whitehouse found a way to generate arbitrarily long 

"cooperative" or totally orthogonal codes for lengths a power 
of 2 [Ref. 3] . 

It will be shown that every complementary code has exactly 
a set of two codes that are totally orthogonal to it. These 
totally orthogonal codes are shown to be transformations of 
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?2 ' + ^ 2 \ 


^ 9^0 + * • • + 


^ ^ -' • 1 - a b 

^ ^ n n 


^2n-l - ^n^l 

Pro. ter. by ter. observation this is recognized to be 
the result of putting a signal series specified by A into a 

frlter .atched to B. c is the cross correlation of A and B 
and can be denoted as follows: 

.T 


C - A0B = diag (a'^B) 

The syrbol x is used to n«an correlation. The abbrevia- 
tron diag is used to denote the row .atrix with the nuntoer of 
ter.s equal to the nuHt,er of diagonals of the .atrix on which 
rt is operating. Each ter. is equal to the su. of the ter.s 
rn the corresponding diagonal. Pro. the previous definition 
of autocorrelation, the same numbers result, but the subscripts 
in the case of autocorrelation range fro. -ntl to n- 1 . i„ the 
case of autocorrelation this presented no difficulty since 
the terms with negative subscripts were equal to the terms 
with positive subscripts. However in the general case it is 
H»re convenient to allow only positive subscripts, particularly 
suice FORTRAN only allows positive subscripts. 

Concern ing Totally Orthoaon.l 
this section the author will develop the theorem 
g that the totally orthogonal codes are transformations 
of the original codes. m the development other theorems 
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regarding linearity, effects of changing the order of cor¬ 
relation, and effects of reversing the codes themselves will 
be proven. 

a. Complementary Correlation 

Let A and B be the two sequences of a complemen¬ 
tary pair, and C and D be another complementary pair. The 
complementary correlation of the first complementary pair 
with the second is defined to be A correlated with B plus C 
correlated with D. This is denoted as: 

= (A® C) + (B® D) 

The symbol is used to indicate complementary correlation. 

Theorem: Complementary correlation is linear. 

That is: 



Where A, B, C, D, E, and F are in the form (a^^, a 2 , • • • , a^) , 
and a and 3 are real constants. Complementary correlation is 
the first sequence of the first pair correlated to the first 
sequence of the second pair, plus the second sequence of the 
first pair correlated to the second sequence of the second 
pair. Thus by the distributive properties of matrices and 
the linearity of the diag operator, the following equations 
can be written: 

(aA + 30 ® E = diag(aA + BC)'^E 

= diag(aA'^E) + diag(Bc'^E) 

= a diag(A'^E) + B diag(c'^E) 
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Similarly: 


(aB + 3D) @ F = a diag(B'^F) + 3 diag(D'^F) 

By adding the above two equations, regrouping, and recognizing 
T T 

that a diag(A E) + a diag(B F) is equal to the complementary 

correlation of a times the pair A, B with the pair E, F; and 
T T 

that 3 diag(C E) + 3 diag(D F) is equal to the complementary 
correlation of 3 times the pair C, D, with the pair E, F; the 
following equation is obtained; 

I“<B> 

Q. E . D. 


b. Theorem: Changing the order of correlation (com¬ 
mutating) results in reversing the correlation. 

To take the transpose of a product of matrices, 
the transpose of each term is taken in reverse order. Thus: 


T 

(B^.A) 


T 


T 

(A^-B) 


The order of the diagonals in the transpose of a matrix is 
reversed. 

A B = diag(A'^*B) 

B0A = diag(B’^.A) = diag (A*^ • B) 


Since the order of the diagonals is reversed, this is the 
same as reversing the cross correlation. 


Q. E. D. 


c. Theorem:, Reversing both codes results in 
reversing the correlation. 


23 


Let 


C 


where A = (a, , a_ , • • • , a. , • • • , a ) 

and B = (b,,b„,•••,b.,•••,b ) 

1' 2' j n 

When correlated the a,b contributes to the first term, and 

i n 

any increase in the a subscript or decrease in the b sub¬ 
script increases the subscript on c. Thus contributes 

to c. and a,b. contributes to c . . Then a.b. contributes 

1 1 3 n+1-] 1 ] 

to • Reversing A and B results the following sequences 


A = (a ,a ,,***,a.,***,a,) 

— n' n-1 ' i' '1 

B = (b ,b T/***,b.,***,b,) 

— n' n-1' D 1 

The i^j^ term in the A sequence is the {n+l-i)^j^ 
term in the A sequence, and the term in the B sequence 

becomes the {n+l-j)^j^ term in the ^ sequence. These terms 
contribute to the correlation terms with the subscript n-i+j. 
Thus the contribution is to the diagonal the same number of 
terms on the other side of the main diagonal corresponding to 
a reverse in the order of terms of the correlation. Since 
this applies to both the first sequences and second sequences 
of complementary correlation, it also follows for complemen¬ 
tary correlation. 


Q. E. D. 


Corollary: Both reversing the correlation and 

changing the order results in the same correlation. This is 
a direct result of applying the preceding two theorems 
together. 
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d. Theorem: A totally orthogonal code can al¬ 
ways be obtained by reversing, exchanging codes, and comple¬ 
menting either one. 

Given: A@B = diag(A *6) 

then by the corollary: 

^ = diag(A'^*B) 


Since the codes are linear: 



B 0 A = 

B 0 A 

= -diag{A'^*B) 

Thus : 

(A0B) + 

(B 0 A) 

= 0 for every term. 

And also 

(A0B) + 

(B 0 A) 

= 0 for every term. 


It should be noted that the above result is 
not dependent upon A and B being a complementary pair. The 
cross correlation will be zero as long as both sequences are 
the same length. Even if each term in the sequence is real 
or even complex it will be true. 

2. Examples 

A few examples are given below. Finding the ortho¬ 
gonal codes of a given code of length four is done in a step 
by step procedure to demonstrate the techniques. The more 
general case of finding the totally orthogonal complementary 
pair sequence pair of length n is done and the complementary 
correlation is shown to be zero. 

All complementary codes of length four are grouped 
into eight groups of a code and its complement and the two 
codes totally orthogonal to them. For any general code there 
are 64 transformations. These transformations are shown to 
partition the 64 transformations into 16 groups, each 
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containing a code, its complement, and the two codes to¬ 
tally orthogonal to them. 

a. Complementary Pair of Length Four 

As the first example, a complementary pair se¬ 
quence of length four whose first sequence is A = +1 +1 +1 -1 
and whose second sequence is B = +1 +1 -1 +1 is considered. 

To get the totally orthogonal sequences using the preceding 
theorem, this pair must be reversed, exchanged and either se¬ 
quence complemented. Underlining is used to denote the 
reverse; while a line over the symbol denotes the complement. 
Reversing: 



A 

= 

-1 

+1 

+1 

+1 


B 

= 

+ 1 

-1 

+1 

+1 

Exchanging: 








B 

= 

+1 

-1 

+1 

+1 


A 

= 

-1 

+ 1 

+1 

+1 

Complementing 

either 

one: 



B 

= 

-1 

+1 

-1 

-1 


A 

= 

-1 

+1 

+1 

+1 

Or: 








B 

= 

+ 1 

-1 

+1 

+1 


A 


+1 

-1 

-1 

-1 


The first one will be demonstrated to be totally 
orthogonal. That is: 

A B 



From the definition of complementary correlation and the 
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representation of correlation presented previously, the fol¬ 
lowing equation can be written; 


A B T- T 

= A@B + B@A = diag(A B) + diag(B A) 
B A “ 


T— 

A B = 


+1 
+ 1 
+ 1 
-1 


(-1 +1 -1 - 1 ) = 


.T;t, 


-1 +1 -1 -1 
-1 +1 -1 -1 
-1 +1 -1 -1 
+1 -1 +1 +1 


A®B = diag(A B) = (-1 -2 -1 0 +1 -2 +1) 



‘+1' 


~-l 

+ 1 

+ 1 

+ 1~ 


+1 


-1 

+1 

+ 1 

+1 

B^A = 

-1 

(-1 +1 +1 +1) = 

+ 1 

-1 

-1 

-1 


+ 1 


-1 

+ 1 

+1 

+1 


B0A = diag(B A) = (+1 +2 +1 0 -1 +2 -1) 

Thus: 


A B 

^c 

B A 


= 0 


Since complementary correlation is a linear operation 
multiplication by a minus one can be done. 

Then: 




B 

A 


0 


Thus the second is also shown to be totally orthogonal to the 
original code. 

b. Complementary Pair of Length n 

As an example the above theorem will be applied 
to a complementary pair binary coded sequence length n where 
each term is either a plus or minus one. Let A and B 
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designate the first and second sequence of the complementary 
pair. 

A@A + B(^B = (0,***,0,2n,0,***,0) 

The totally orthogonal codes in this case are 
B, A and B, A. This first pair of these when written out 
term by term is: 


B = (-b , -b , , * • * , -b,) 

— n' n-1 1 

A = (a ,a T , • • • ,a,) 

— n n-1 ' 1 

This first pair of these is shown to be totally orthogonal 
to the original pair, but since complementary correlation 
is linear to show that the second is totally orthogonal can 
be very easily done by observing the second pair is just 
minus one times the first pair as in the previous example. 


T — 
A ‘B 


B 


n 




L^J 


= diag(A'^B) + diag(B'^A) 


(-b , -b , , • •t,-b,) = 
n' n-1 1 


-a, b 

-a. b 

,•••-a,b,! 

1 n 

1 n 

-1 1 1| 

-a^b 

-a„b 

, • • ‘-a^b,. 

2 n 

2 n 

-1 2 1' 

• 

• 

• i 

• 

• 

• 1 

• 

• 

• 

-a b 

-a b 

•••-a b ‘ 

n n 

n n 

-1 n Ij 


B X A = diag(A B) 


T 

B ‘A 


n 


(a / a , , 
n' n-1 




b, a 

b, a 

• ♦ ♦b,a 

1 n 

1 n 

-1 1 , 

b_ a 

b_a 

- • •‘b^a 

2 n 

2 n 

-1 2 , 

• 

• 

• 

• 

• 

• 

• 

• 

• 

b a 

b a 

,••"b a 

n n 

n n 

-1 n 
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# I 


Thus when the diagonals from above matrices are 
added together, for every minus number in the upper diagonal 
the positive term is in the same diagonal in the second ma¬ 
trix. The complementary correlation, which is the sum of the 
two individual cross correlations, is zero for every term. 

c. Codes of Length Four 

An exhaustive computer search was done to find 
all totally orthogonal codes of length four. The thirty-two 
possible complementary pair sequences including all possible 
transformations were read into the computer. Only the auto- 
correlation, the totally orthogonal codes, their correlation, 
and their complementary correlation were printed out. See 
the results in Table I. 

The results can be summarized by saying for each 
code and its complement a totally orthogonal code and its 
complement were found. The thirty-two possible codes includ¬ 
ing all transformations could be partitioned into eight 
groups of four. It can be easily verified that each to¬ 
tally orthogonic to another is the reverse, interchange, 
complement of one. 

d. Transformations of Complementary Codes at Any 

Length. 

Previously Golay described 64 transformations that 
were possible for each complementary code. Jauregui [Ref. 2] 
showed that they were in fact a closed group, identified all 
possible operators, and demonstrated an operation multiplica¬ 
tion table. Defining operations as in Table II, and using the 
Jauregui operation multiplication table it can be 
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TABLE I 


EXHAUSTIVE LIST OF TOTALLY ORTHOGONAL 
COMPLEMENTARY PAIRS OF LENGTH FOUR 


Standard 

Code 


Complement 

Code 


Totally 

Orthogonal 

Code 


Complement 

Totally 

Orthogonal 

Code 


+1 

-1 

-1 

-1 

-1 

+1 

+1 

+1 

+1 

+1 

-1 

+1 

-1 

-1 

+1 

-1 

-1 

+1 

-1 

-1 

+1 

-1 

+1 

+1 

-1 

-1 

-1 

+1 

+1 

+1 

+1 

-1 

+1 

-1 

-1 

-1 

-1 

+1 

+1 

+1 

+1 

-1 

+1 

+1 

-1 

+1 

-1 

-1 

-1 

-1 

+1 

-1 

+1 

+1 

-1 

+1 

-1 

-1 

-1 

+1 

+1 

+1 

+1 

-1 

-1 

-1 

-1 

+1 

+1 

+1 

+1 

-1 

+1 

+1 

-1 

+1 

-1 

-1 

+1 

-1 

-1 

+1 

-1 

-1 

+1 

-1 

+1 

+1 

+1 

-1 

-1 

-1 

-1 

+1 

+1 

+1 

-1 

-1 

-1 

+1 

+1 

+1 

+1 

-1 

+1 

-1 

+1 

+1 

-1 

+1 

-1 

-1 

-1 

-1 

+1 

-1 

+1 

+1 

-1 

+1 

+1 

-1 

-1 

-1 

-1 

+1 

+1 

+1 

-1 

+1 

-1 

-1 

+1 

-1 

+1 

+1 

+1 

+1 

+1 

-1 

-1 

-1 

-1 

+1 

+1 

-1 

-1 

-1 

-1 

+1 

+1 

+1 

-1 

-1 

+1 

-1 

+1 

+1 

-1 

+1 

-1 

-1 

+1 

-1 

+1 

+1 

-1 

+1 

+1 

+1 

+1 

-1 

-1 

-1 

-1 

+1 

+1 

-1 

-1 

-1 

-1 

+1 

+1 

+1 

-1 

+1 

-1 

-1 

+1 

-1 

+1 

+1 

-1 

+1 

-1 

-1 

+1 

-1 

+1 

+1 

-1 

+1 

+1 

+1 

+1 

-1 

-1 

-1 

-1 

-1 

-1 

+1 

+1 

+1 

+1 

-1 

-1 

-1 

-1 

-1 

+1 

+1 

-1 

+1 

-1 

-1 

+1 

-1 

+1 

+1 

-1 

+1 

-1 

+1 

+1 

+1 

+1 

-1 

-1 

-1 

-1 

-1 

-1 

+1 

+1 

+1 

+1 

-1 

-1 

+1 

-1 

-1 

+1 

-1 

+1 

+1 


Note: Every complementary code of length four appears in the 


table once. In its row is contained its complement and the 
two complementary codes that are totally orthogonal to it. 
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TABLE II 


ELEMENTS OF THE UNORDERED OPERATIONS GROUP 


Sequence A=a, a„a,‘**a -.a • 

^ 123 n-1 n 




II=a a „ 
n n-2 


n-1 


Sequence 




n-1 n 


III=b^b3b5 


IV=b b „ 
n n-2 


n-1 


I 

= 

(I II 

III 

IV) 

R 

z= 

(I II 

IV III) 

^1 

= 

(1 II 

III 

IV) 

Q 

= 

(II I 

IV III) 

^2 

= 

(I TT 

III 

iv) 

p 

= 

(iT I 

iv III) 

^1 

= 

(II I 

III 

IV) 

0 


(iT I 

III IV) 

^2 


(I II 

IV 

III) 

N 

— 

(II I 

III IV) 

T 

= 

(II I 

IV 

III) 

M 

— 

(I rr 

IV III) 

^1 


(I II 

III 

IV) 

L 


{TT T 

III IV) 

^2 


(I II 

III 

IV) 

K 


{T TT 

IV III) 

C 

= 

(I iT 

III 

IV) 

J 

= 

(iT I 

IV III) 

z 

= 

(II I 

III 

IV) 

H 

=r 

(I II 

IV III) 

Y 

= 

(iT I 

ill 

IV) 

G 

= 

(iT r 

IV III) 

X 

= 

{TT I 

III 

iv) 

F 


(I iT 

III IV) 

W 

= 

(II I 

III 

rv) 

D 


(T II 

III rv) 

V 

= 

(I II 

IV 

III) 

0 

= 

(II I 

IV III) 

U 

= 

(I II 

rv 

III) 

B 

= 

(iT I 

IV III) 

S 


(I TT 

iv 

E 

III) 

ir 

(III IV I 

II) 

(II T 

iv III) 


Note: Combining these operations with the E (exchange) 

operation provides the additional 32 transformations 
for a total of the 64 possible transformations. 
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TABLE III 


GROUPS OF OPERATIONS PRODUCING TOTALLY ORTHOGONAL PAIRS 


Operation 

Complement 

Operation 

Totally 

Orthogonal 

Operation 

Complement 

Totally 

Orthogonal 

Operation 

Operations 
Performed left 
Hand Member 

E C^ 

I 

C 

EG 

E0 

0 

0 

0 

0 

^2 

^1 

EB 

EH 

0 

0 

0 

1 

^1 

^2 

ET 

EJ 

0 

0 

1 

0 

D 

F 

EQ 

EP 

0 

0 

0 

1 

1 

^1 

L 

EM 

EH 

0 

1 

0 

0 

X 

Z 

ER 

EU 

0 

1 

0 

1 

0 

N 

ET2 

EK 

0 

1 

1 

0 

w 

Y 

EV 

ES 

0 

1 

1 

1 

El 

EC 

. G 

■0 

1 

0 

0 

0 

EA2 

EA^ 

B 

71 

1 

0 

0 

1 

ECj^ 

EC2 

T 

J 

1 

0 

1 

0 

ED 

EF 

Q 

P 

1 

0 

1 

1 

ET^ 

EL 

M 

H 

1 

1 

0 

0 

EX 

EZ 

R 

U 

1 

1 

0 

1 

EO 

EN 

T 

2 

K 

1 

1 

1 

0 

EW 

EY 

V 

S 

1 

1 

1 

1 


Note: Each operator appears only once and there is another 

operation which produces a code which is the complement to 
the code produced by the given operation. There are two 
operations that produce codes totally orthogonal to the code 
produced by the first operation and its complement. 
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demonstrated that the 64 transformations are partitioned into 
16 subgroups of code, complement, orthogonal code, and 
complement to the orthogonal code. See Table III. 


3. Linear Homogeneous Equation Approach 

Another way of expressing correlation is by a matrix 
multiplication as shown below: 


Let 


C = c^,C2,*--,c^ 


A = 


n 


TP — -p -p ••••p 

!■ '^2n-l 


C 09 A = F 

This correlation can be represented as shown below: 


(Ci,C2, 




a a - ■ 
n n-1 

0 a a , 
n n-1 


a^ 0 • • • 0 
a^ 0 


0*r*0 a a 

n n-1 


(fl,f2,••*/f2n-l^ 


m 


= F 


The subscript m is used to denote a matrix which is more than 
just a row matrix. This notation is used for convenience to 
show the form of the matrix equations. 

For complementary correlation, let the n-dimensional 
vectors C and D represent one code and the n-dimensional 
vectors A and B represent the second code. Let channel one 
correlation be the vector F of dimension 2n-l and channel two 
correlation similarly be the vector G. Let the vector H be 
the complementary correlation or the sum of F and G. 
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c 


I—1 


u 


o 


x 

•H rH 

M I 

-P d 
fd cN 



AJ U} 

o 

fd P • 
^ U) 

-P d d 

<N fc 

Q) d 

-P U) r-H 

o fd o 
^ s: o 
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For convenience the transpose of both sides is taken. 



To find the totally orthogonal complementary pair 
involves solving the linear homogeneous equation with 2n 
unknowns and 2n-l equations- By homogeneous is meant that the 
H vector is the null vector. This can be represented by the 
matrix equation: 


A X = 0 

m 

A is the matrix with 2n-l rows and 2n columns. The answer 
m 

to this problem has been worked in many books [Ref. 4]. 


X. = k A 
1 ' m. ' 

1 

The matrix A is (-1) matrix with the i,, column 

m. m th 

1 

crossed out, and k is an arbitrary constant, 
a. Example: 

Using the linear homogeneous equation approach, 
find the complementary pair totally orthogonal to the comple¬ 
mentary pair, A = (1,1), B = (1,-1). 
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— 




r *1 

+1 0-1 0 


^1 


0 

+1 +1 +1 -1 


^2 

II 

0 

0 +1 0 4-1 


<^1 


O 


1 CN 

_1 











0 

-1 

0 





-41 

-41 

-1 





-41 

0 

-41 





-41 

-1 

0 






-41 

-41 

-1 



2 

^1 



0 

0 

+ 1 




^2 

= k 





= k 

-2 

d. 



-41 

0 

0 




1 

d. 



-41 

-41 

-1 



2 

2 



0 

-41 

-41 



2 



-41 

0 

-1 





-41 

-41 

-41 





0 

-41 

0 










Restricting c^^, C 2 f 

di. 

and d^ to -41, 

and 

-1 the possible 

solutions are: 





First solution 

(Cl 

^2 ^1 ^2^ “ 

( + 1 

-1 -41 -41) 

Second solution 

(Cl 

^2 ^1 ^2^ 

(-1 

-41 -1 -1) 


Clearly the linear homogeneous equation can 
become very difficult to solve because of the size of the 
determinants even for comparably short codes. For the case 
where n=4 the solution involves solving eight seven dimen¬ 
sional determinants. However no matter how we get a solution, 
a constant times the solution is also a solution. Thus re¬ 
stricting ourselves to +1 and -I's the solution found previ¬ 
ously (the reverse, exchange, and the complement of one) and 
its negative are both solutions. 
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The question arises as to whether it is possible 
that there could be another solution orthogonal to both the 
original complementary pair and to the solution that was 
found. 

b. Theorem: 

There is no complementary code totally orthogonal 

to both the original code and to one orthogonal code already 

found by reversing, exchanging, and complementing either one. 

The correlation of both the original code and 

the orthogonal code found with a new code can be hypothesized. 

Let the solution be (c,,***,c ,d,,***,d ). Then being ortho- 

1 n l n ^ 

gonal to the first involves 2n-l equations and being ortho¬ 
gonal to the second involves another 2n-l equations. These 
can be gathered together into one matrix. 
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n 


0 

b. 


n 


0 

a 


n 


0 a. 

0 0 

0 0 


0 

0 

b. 


n 


0 

0 


0 

0 


n 


0 

0 

0 

0 

b. 


n 


0 

0 

0 


n 


0 

b. 


n 


0 

0 

0 

0 

a 


n 


0 a- 

0 0 

0 0 


0 

0 

b. 


n 


0 

0 

0 

0 

a. 


n 


0 

b 


n 


0 0 
0 0 


0 

0 

b 


n 


0 0 
0 0 


n 


0 

0 

0 

b 


n 


0 


n 


0 

0 

0 

0 

b 


n 


0 0 0 0 

a^ 0 0 0 

32 0 0 

* ^2 ^1 ° 

a2 3] 

a • • a, 

n ^ 


n 


n 


This is in the form of A X = 0. In this case 

m 

there are 2n unknowns and 2(2n-l) equations. has rank 

less than or equal to 2n, because it only has 2n columns. 

remains to show A has rank 2n. 

m 

From linear algebra the rank of the product of 
matrices is not greater than any of its factors [Ref. 5]. 
For example: 


A ‘B 
m m 


= C 


m 


rank(C ) < min (rank(A ), rank(B )) 

m — m ' m 


It 
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Thus a corollary of this is that the rank of the product of 
a matrix and its transpose is not greater than that of the 
matrix itself. 

T 

rank(A • A ) < rank(A ) 


T 

A .A = 
m m 


n 


^ 0 


0 

b 


' 0 a 


n 


m m' 

0 - 

0 • 

* 

0 . 

a^^O • * 


b 0 
n 


0 b. 


0 

b 


n 


n 


0*•*0 a ^ 
• • ; 0 
’ • 0 ; 


a 0 
n 


Ob 

n 


b^O 


0 a. 


n 


a 0 • • 

• 0 b 0 • 

• ‘0 

n 

. n 

• 

- 

• 

. 0 ; 

• 

' . 0 

a, . . 

- a b, • 

• ■ b 

1 

n 1 

n 

0 . 

■ 0 • 

• 

• 

• 


0* • *0 

a, 0 • • • 

0 b. 


1 

1 

bj^O • ' 

'0 a^^O • 

• -0 
• 

* * 

• • • 

• 


. 0 • 

* 0 

b . . • 

b, a . 


n 

1 n 

1 

0 . 

• 0 . 

. 


0***0 b a 

n n__ 

When multiplication is done to get the main diagonal terms 
of the product it is just the autocorrelation of a complemen¬ 
tary pair sequence. Any other position is either a shifted 
correlation of a complementary pair or correlation with the 
reversed, exchange, complement of one which is totally ortho¬ 
gonal. In either of the last two cases the result is zero, 
but the main diagonal is 2n in each position. 

Therefore: 

a'^-A = 2n 

mm 2n 

T 

The rank of A -A is 2n. The rank of A is greater or equal 
mm m 

to 2n, but we previously showed that it was less or equal to 


2n. Therefore the rank of A is 2n. 

m 
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V 


Since the rank of A is 2n, there is no non- 

m 

trivial solution to the equation, and no code totally ortho¬ 
gonal to both the original code and one of the totally 
orthogonal codes can exist. 

Q. E. D. 
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III. PROPOSED COMMUNICATIONS SYSTEMS 


A matched filter for a given signal can be built either 
by analog or digital means to a complementary coded signal. 
The signal itself and its complement can be used as the two 
signals mark and space, or one and zero. A number of varia¬ 
tions are proposed, evaluated for different applications. 

The first category is variable length coding system dependent 
upon noise, and the second system proposed is a multiplexing 
system based on complementary coding. 

A. VARIABLE LENGTH CODING 

1. Two Length Codes for Same Filter 

Dr. S. Jauregui, Jr., mentioned the possibility of 
having a communication system based upon using a long comple¬ 
mentary pair filter for high noise levels and when noise 
levels were lower a shorter pulse sequence could be used, but 
the filter would remain the same. 

Example: 

Let A = +1 +1 

B = +1 -1 

AB = +1 +1 +1 -1 

AB = +1 +1 -1 +1 

Time reversing the bottom code, 

+1 +1 +1 -1 

+ 1 -1 +1 +1 
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This is a standard technique for generating longer codes, and 
the time reversing of the second code is a standard trans¬ 
formation, so this is a complementary sequence. 

Autocorrelation 

+1 +1 +1 -1 +1 + 1 + 1-1 -10141 0-1 

+1 -1 +1 +1 +1 -1 +1 +1 1 0-1 4-1 0 1 

Complementary correlation 0008000 

A shorter sequence in the same filter is: 

+1 +1 +1 +1 +1 -1 -1 0 2 2 1 

+ 1 -1 +1 -1 +1 +1 1 0-2 2-1 

Complementary correlation 00040 

Dr. Jauregui demonstrated the same phenomenon by hand for 
length of 8. Investigation of the following codes was con¬ 
ducted using the IBM 360: 

Length 16 with code lengths 8,4,2 
Length 40 with code lengths 20,10 
Length 104 with code lengths 52,26 
The method of generating the codes is: 


B. 


AfEi 

Vl 



^ 2^2 

A2B2 



etc. , 


where B is the complement of B. At the final step the lower 
code is time reversed. 

The digital computer flow chart is in Fig. 6. The re¬ 
sults are summarized below. In each case n is the signal 
code length. 
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*Dashed line is not normally used. Exit is normally made 
from the middle loop. 


Figure 6. Variable Length Pulse Digital Computer Flow Chart 
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Half filter code length correlation 


Mark 0 ^, 

Space 


■ ■ ■ ''’Sn-l 


Quarter filter 

length correlation 


Mark ^i'***'°n-l' 



• • • 0 

'^5n-l 

Space 0i'“*/0n-l' 

■"'“n+l'' 


• • • n 

'^5n-l 

One-eighth filter length 

correlation 


All terms are zero 

except as 

noted below: 


Mark 

Term Value 

Space 

Value 


n 

-n 

+n 


2n 

+n 

-n 


3n 

+n 

-n 


4n 

+n 

-n 



For codes of length one-half of the filter length, a 
single output pulse is generated. For shorter length codes 
the multiple pulse output makes the decoding somewhat ambigu¬ 
ous if the signals are sent sequentially immediately following 
the previous one. 

An explanation of why this phenomenon v/orks so well for 
codes of half length is given below. Taking the final two 
steps of building the code (underlining means reversing): 

Given 


A 

Generating 

A B 

Reversing 

A B 

Simplifying 

A B 

B 

a new code 

A B 

second co(?e 

A B 

- ^ 

B A 

- - 
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I 


1 

1 




# 

■* 








The first part of the code is totally orthogonal to the 
second part. This can be demonstrated by showing that the 
second part is the exchange, reverse, complement of one 
transformation. 


m 4 

A 

Exchange 

ft i 

B 

Time inverse 

• m 

B 

Complement one 

• 4 

B 

B 

• « 


A 

V • 


A 

■ « 


A 

• « 


Thus the first part of the code is orthogonal to the second 

c c * 

part. Consider the code . Let (p,) be its totally 

orthogonal code. Since correlation is linear the following 
can be written: 


m m 





^ m 









c 

®c 

C 

C ' 

= 

c 


C 

0 

+ 

c 

®c 

0 

c 

D 


D 

D'_ 




D 

0_ 


D 


0 

D' 


Thus: 


C 

®c 

c c' 

D_ 


D D' 



, 0 T ,n, 0 . T , 

n-1 n+1' 


* '°3n+l] 


2. Totally Orthogonal Code Scheme 

By alternately sending the coded signal and the 
orthogonal coded signal, it is possible to achieve maximum 
time separation of the different information signal on the 
two channels. See Fig. 5 for equipment setup. 

The advantage to this system is that if multiple paths 
tend to spread signals over a period of time, the orthogonal 
sets keep the channels separated in time as far as possible. 
For instance in the above case with code lengths of four, the 
first message bit is on the standard code. The next bit on 
this code would come 8 bits later. In the meantime a message 
is decoded on the orthogonal code. 
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The total number of bits of information is as if the 


length of the code was four bits but the separation time on 
each channel is the time to send 8 bits. In the intervening 
time all signal generated clutter is canceled between the two 
channels. In other words, the transmitter can broadcast con¬ 
tinuously, but when each message is reconstructed at the 
received, it is like an impulse without any side lobes which 
in a scatter environment could be falsely interpreted as 
signals. 

B. MULTIPLEXING 

The property that eliminates clutter from the radar sys¬ 
tem allows the improvement in resolution for a series of 
pulses is the orthogonality of the code with a time shift of 
itself. In a communication system this would allow the sender 
to start the next message one bit behind the previous one and 
by superposition of this linear system, it can be seen that 
these bits will be separated at the receiver end. 

For simplicity the synchronous analog ideal receiver will 
be assumed for the uncoded pulse. The white gaussian noise 
will cause gaussian noise at all frequencies. Therefore it 
will have gaussian noise out of the synchronous detector. 

The signal may be put on a carrier, but at the receiver 
a synchronous receiver is hypothesized. For the matched fil¬ 
ter, the impulse response is just the time reverse of the 
signal. Then the normal convolution of the signal and filter 
becomes a correlation of the signal with itself. For con¬ 
venience a square pulse is assumed for the uncoded system. 
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One way to get this impulse response is to use an integrator. 
The two signals are two pulses of equal magnitude and length 
but opposite phase. The initial condition of the integrator 
is zero. Integration proceeds during the length of the 
impulse. At the end of the baud, the value is sampled. The 
integrator is then reset for the next pulse. If it is posi¬ 
tive the decision is made that the pulse was in phase, if 
negative the decision is made that the pulse is out of phase. 

For a longer pulse, say n times the original length, either 
the signal could be sampled n times, the values added together 
and the decision made on the sum or the integrator could be 
sampled at the end of the total length. 

For a coded pulse the signal can be sampled at each part 
and by multiplying either by plus or minus one according to 
the code and summing a decision can be made at the end as to 
which signal was sent. 

Similarly for the complementary code, the scheme for a 
coded pulse is applied to both channels and the results are 
added together. The decision is made based on the sum of 
both channels. 

Since the output of the matched filter itself is gaussian, 
it can be simulated by a gaussian noise generator. Since the 
uncoded matched filter is a part of the complementary code 
matched filter simulation, this can also be done for the multi¬ 
plexed system, and thus provides a valid comparison, between 
the basic uncoded system, the basic complementary system, the 
multiplexed system using a single code, and the full 


47 



multiplexed system. The signal to noise ratio is the power 
ratio measured at the output of the basic matched filter 
system. 

For the uncoded system the noise is assumed gaussian at 
the output of the matched filter. The probability of an 
error is: 

P(error)= P (mark)-P (error mark) + P (space)‘P (error space). 

The probability of error given a mark is the same as the 
probability of an error given a space, and the probability 
of a mark plus the probability of a space is one. Then: 

P (error) = 1 - F(x) 

Where F(x) is the cumulative distribution function of the 
normal or gaussian probability, function x is the square root 
of the signal to noise power ratio. 

A complementary code of length sixteen was used to en¬ 
code a message of length one thousand bits and the number of 
errors were recorded. The signal to noise ratio required for 
any error percentage was 32 times lower, but two channels 
were needed and the flow of information was 16 times slower 
(See Fig. 7). 

Full multiplexing using both the standard and totally 
orthogonal codes was conducted for codes of length 128, 64, 
32, 16, 8, 4, and 2. All produced results in close agreement 
with the basic uncoded and the theoretical curve based on 
gaussian distribution for the uncoded system (See Fig. 8). 

If higher data rate is desired the series of pulses in 
the code can be compressed into a shorter time. This results 
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in a wider bandwidth, but using the horizontal multiplexing 
technique, the bandwidth can be fully utilized. 

Depending upon the noise background the amount of signals 
multiplexed onto a complementary channel can be limited to get 
any desired accuracy. If very low noise level is present 
another method may be used to again increase the transmitted 
information. 

An additional fifty percent of information can be trans¬ 
mitted over the bandwidth and channel if the noise level is 
sufficiently low. Three bits are taken at a time either from 
three sources or from one source at a higher rate. This is 
converted to a base three number and mapped into the plus one, 
minus one, or zero and encoded. At the receiver it is re¬ 
ceived as plus one, minus one, or zero on both standard and 
the totally orthogonal filter, mapped into a ternary system 
converted to binary. 


bits 

Ternary 

Standard 

Filter 

Orthogonal 
FiIter 

Ternary 

Decoded 

000 

00 

-1 

-1 

00 

000 

001 

01 

-1 

1 

01 

001 

010 

02 

-1 

0 

02 

010 

Oil 

10 

1 

-1 

10 

Oil 

100 

11 

1 

1 

11 

100 

101 

12 
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0 

12 

101 

110 

20 
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20 

110 

111 
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21 

111 

> signal 22 

0 
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22 

No signal 


Mapping 

Reverse Mapping 
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T-423 is a digital to analog trunk line. By calling the 
subroutine DAC a voltage specified from the digital computer 
is set on that trunk line. T-500 is an analog to digital 
trunk line. By calling the subroutine ADK the analog signal 
present at this trunk line is sampled, converted to digital, 
and stored in the digital computer under the variable specified. 

Integrator A-OOl is used as a track-hold network; holding 
the voltage while in compute mode, and tracking the voltage 
when the analog computer is in reset. The capacitor selected 
for A-OOl was the smallest available, .001 microfarads to 
allow very rapid reset times. 

The noise generator was fed into an amplifier to maintain 
isolation under all settings of P-000. 


Figure 9. 


Analog Computer Diagram 




























Figure 10. Digital Computer Flow Chart of Uncoded Pulses 
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Figure 11. Digital Computer Flow Chart of Unmultiplexed 
Complementary Coded System 
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IV. CONCLUSION 


1. Totally orthogonal complementary binary coded sequences 
were found to be the exchange, reverse, complement one 
transformation of the basic code. The basic code and its 
complement were both proven to be totally orthogonal to 
the orthogonal code and its complement. 

2. These totally orthogonal codes were developed with an eye 
toward communications, however applications also exist in 
sonar, and possibly radar. If nothing else it would be 
possible to double the unambiguous range for any given 
repetition frequency by alternately sending the standard 
code and the totally orthogonal code and receiving both 
on the standard code and the orthogonal code filters. 

3. Depending on the noise level the code length can be ad¬ 
justed as necessary to allow the transmitter to send how¬ 
ever much energy per bit necessary to communicate reliably. 
However the data rate has to go down during noisy times 

to maintain the necessary fidelity of information. 

Virtually any rate of multiplexing can be done on the same 
bandwidth to obtain the error probability required. 

4. In the system simulated white gaussian noise was assumed, 
but various burst noise levels could be simulated, and if 
the duration of the burst were short, the signal can still 
be decoded correctly. Since the signals are really spread 
over a long time on the two channels by deleting a few 
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bits here and there would not appreciably affect the net 
result. This method could be very advantageous if used 
properly. 

5. Advantages of this multiplex system would be the ease 

of splitting the signal up using digital techniques, and 
without the need for guard bands as in the frequency 
multiplex system. An advantage of this system over the 
time multiplex system is the bit error correction in a 
burst channel noise situation. Among the inefficiencies 
of needing several levels of signals in the output trans¬ 
mitter. The need for two channels might be overcome by 
sending the two signals in quadrature, or even time multi¬ 
plexing if the signal is to be digitally processed. 
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APPENDIX A 


SIMULATED COMMUNICATION SYSTEM PROGRAMS 

Three digital programs are included in this thesis. The 
first two correspond to the flow charts in Figs. 10 and 11 
for the vincoded system and simple coded systems. The third 
is the full multiplexed system using both the totally ortho¬ 
gonal code and the one bit delay. 

All programs were debugged on short codes and the results 
of each bit was recorded, but after debugging only the sum¬ 
mary of each rxon was printed. This was accomplished by using 
the X in the comment column of the write cards which were to 
be eliminated. By using the control card "FORTRAN LS,GO,X" 
these are compiled, but when "FORTRAN LS,GO" is used, these 
cards are treated as comments. 

The length of code can be changed by changing input data 
and three cards: (1) DIMENSION statement, (2) DATA statement 
for KDIM, and (3) 8501 FORMAT statement. Similarly the mesr- 

sage length can be changed by changing the (1) DIMENSION 
statement, (2) DATA statement for MESDIM, and (3) 8500 FORMAT 

statement. 

On the third program it was necessary to change some 
simple variables into arrays. The message of the full multi¬ 
plexed system requires one bit less than the sum of the 
itessage length and the code length. On fairly short messages 
and long codes this could distort the amount of noise power 
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in the calculation. This program corrected for this to make 
the proportion the same as for the case that the message is 
much longer than the code length. 
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PRINCIPAL VARIABLES USED 


CODE 

compleimentary code(s) in real numbers. 

DBSN 

signal to noise level in decibels. 

lERR 

increased by one if there is an error and is zero 
otherwise. 

KDIM 

dimension of complementary code. 

KODE 

complementary code array in integers. 

KPOT 

array of pot settings for various runs. 

KPOTDM 

dimension of KPOT. 

MES 

message composed of O's and I's read in for con¬ 
venience and then changed to -I's and +l's. 

MESDIM 

length of the message 

NERR 

number of errors. 

NUMBIT 

used in totally orthogonal multiplexing due to two 
channels of message and is equal to 2 times riESDIM 

OMES 

plus or minus one based on received signal. 

QMES 

array of preceding message bits, used in multiplex 
encoding. 

RMES 

received signal. 

RMEST 

s\im of received coded bits, on which the final 
decision of whether the reception is +1 or -1. 

SIGXN 

signal to noise ratio of each bit. 

SIGXNT 

signal to noise ratio of total received message. 

SMES 

message sent (real numbers). 

SPWR 

signal power each bit (normalized to 1 ohm). 

SPWRT 

signal power in the entire message. 

XNOISE 

the difference between SMES and RMES, the noise 
voltage for each bit. 
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XPWR 

XPWRT 


noise power, i.e., XNOISE squared, 
noise power for the total message. 
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WRITE(6,9000) MES(I),OMES,SIGXN 
IF (lERR.EQ.O) GO TO 200 
WRITE(6,9001)NERR 
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